The evolution from BCS to BEC superfluidity in the presence of disorder 



O 
(N 

o 

Q 
o 



C/2 



o 



X 



Li Han and C. A. R. Sa de Melo 
School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA 

We describe the effects of disorder on the critical temperature of s-wave superfluids from the BCS 
to the BEC regime, with direct application to ultracold fermions. We use the functional integral 
method and the replica technique to study Gaussian correlated disorder due to impurities, and we 
discuss how this system can be generated experimentally. In the absence of disorder, the BCS regime 
is characterized by pair breaking and phase coherence temperature scales which are essentially the 
same allowing strong correlations between the amplitude and phase of the order parameter for super- 
fluidity. As non-pair breaking disorder is introduced the largely overlapping Cooper pairs conspire 
to maintain phase coherence such that the critical temperature remains essentially unchanged, and 
Anderson's theorem is satisfied. However in the BEC regime the pair breaking and phase coherence 
temperature scales are very different such that non-pair breaking disorder can affect dramatically 
phase coherence, and thus the critical temperature, without the requirement of breaking tightly- 
bound fermion pairs simultaneously. In this case, Anderson's theorem does not apply, and the 
critical temperature can be more easily reduced in comparison to the BCS limit. Lastly, we find 
that the superfluid is more robust against disorder in the intermediate region near unitarity between 
the two regimes. 

PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Fk 



I. INTRODUCTION 

Ultracold atoms are special systems for studying super- 
fluid phases of fermions or bosons at very low tempera- 
tures, because of their unprecedented tunability. In par- 
ticular, ultracold fermions with tunable interactions were 
used to study experimentally the so-called BCS-to-BEC 
evolution, where fermion superfluids were investigated 
as a function of the interaction parameter (scattering 
length) . The study of superfluidity in ultracold fermions 
has additional promising research directions which in- 
clude the BCS-to-BEC evolution in optical lattices 
and the effects of disorder during the BCS-to-BEC evo- 
lution, which would allow the very important study of 
the simultaneous effects of interactions and disorder at 
zero 0, Q and finite temperatures @, 0| • 

In ordinary condensed matter (CM) systems the con- 
trol of interactions is not possible, and the control of dis- 
order caused by impurities is very limited because their 
concentrations can not be changed at the turn of a knob. 
However, in ultracold atoms it may be possible to create 
disordered impurity potentials, where control over the 
impurity concentration and potential could occur in con- 
junction with the control over atom-atom interactions. 
Impurity type potentials have been very recently created 
experimentally in the context of thermometry for spin- 
dependent optical lattices [sj. However, we envision the 
creation of another type of impurity potential where there 
are two types of atoms, e.g., a mixture of atoms with un- 
equal masses The first type of atom (the heavier 
one) is tightly trapped by an optical lattice in a low fill- 
ing fraction configuration, such that the locations of the 
atoms for any given realization is unknown (random). 
The second type of atom (the lighter one) is trapped by 
a harmonic potential, but does not feel the optical lattice 
potential of the first type which exists in the same spa- 



tial region. In this case, the lighter atoms only interact 
with randomly distributed heavier atoms which are then 
seen as randomly distributed impurities, and both the 
density of impurities (heavier atoms) and the scattering 
potential between the impurities (heavier atoms) and the 
lighter atoms can be controlled. This kind of situation is 
illustrated schematically in Fig. [1] 




FIG. 1: Schematic illustration of the heavier atom (e.g. *°K) 
impurities trapped in an optical lattice and the trapped cloud 
of lighter atoms (e.g. ®Li) which become superfluid at low 
temperatures and low impurity concentration. 



Another way of introducing disorder in ultracold atoms 
relies on the use of laser speckles or lasers with incom- 
mensurate wavelengths, which were used to study the 
phenomenon of Anderson localization in ultracold Bose 
atoms im ■ The disorder introduced in this way is 
very different from impurity disorder in two important 
aspects. First, laser speckles or lasers with incommensu- 
rate wavelengths produce a disorder potential landscape 
that have many ups and downs [valleys and mountains) 
over a characteristic length scale comparable to the laser 
wavelength, while in the case of impurity disorder consid- 
ered here, the disorder potential is purely repulsive, that 
is, has only ups {mountains) located at random positions. 
Second, in laser speckles, the disorder is exponentially 
correlated, while in the impurity problem the disorder is 
Gaussian correlated. These two major differences may 
lead to results that are both quantitatively and qualita- 
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tively different, enforcing the idea that not ah types of 
disorder produce the same results. Two interesting re- 
view articles have emerged recently covering mostly the 
effects of disorder in ultracold Bose atoms O [l^ ■ 

In this manuscript, we describe the critical tempera- 
ture of three dimensional (3D) s-wave Fermi superfluids 
from the BCS to the BEC limit as a function of disorder, 
which is independent of the hyperfine states of the atoms 
and is created by randomly distributed impurities. How- 
ever, the effects of disorder during the BCS-to-BEC evo- 
lution in higher angular pairing 14] or for spin-dependent 
(hyperfine-state-dependent) impurity potentials [6] can 
also be explored. 

Our main results are as follows. First, in the BCS 
limit the amplitude and the phase of the order parame- 
ter are strongly coupled, such that pair breaking and loss 
of phase coherence occur simultaneously. In this case, 
the critical temperature is essentially unaffected by weak 
disorder, since the disorder potential is not pair-breaking 
and phase coherence is not easily destroyed in accordance 
with Anderson's theorem [l^. Second, in the BEC limit 
the breaking of local pairs and the loss of phase coher- 
ence occur at very different temperature scales. In this 
case, the critical temperature is strongly affected by weak 
disorder (in comparison to the BCS regime), since phase 
coherence is more easily destroyed without the need to 
break local pairs simultaneously, and Anderson's theo- 
rem does not apply. Third, we find that superfluidity is 
more robust to disorder in the intermediate region near 
unitarity between the BCS and BEC regimes. 

The rest of the manuscript is organized as follows. In 
section |lTl we discuss the Hamiltonian for the impurity 
problem and the corresponding impurity potential. In 
section IIIIl we describe the effective action in the pres- 
ence of disorder near the critical temperature, where the 
order parameter for superfiuidity is small. In section HVl 
we apply the replica trick to obtain the critical tempera- 
ture as a function of disorder and interaction parameters. 
We calculate the thermodynamic potential first and then 
the number equation in the presence of disorder. In sec- 
tion |Vl we present analytical results for weak disorder 
both in the BCS and BEC limits. In section IVTl we ana- 
lyze the mean free paths in limiting cases and relate them 
to the parameter describing disorder. Before concluding, 
we compare qualitatively speckle and impurity potentials 
in section fVIII Finally, we summarize our conclusions in 
section IVIIIl 



II. HAMILTONIAN 

In order to study the effects of disorder in Fermi su- 
perfluids from the BCS to the BEC regime we start with 
the real space Hamiltonian density for three dimensional 



s-wave superfluids (set h = 1) 

H(x) + ^dis(x)') ^.(x) + ;7(x), 

(1) 

where il'ti^) represents the creation of fermions with 
mass m and hyperfine state (spin) a at x, Vdis(x) is the 
disorder potential, and /i is the chemical potential. In 
addition, the interaction term 

U{^) - + y dx'y(x,x')V4(x')V'I(x)^i(x)V't(x') (2) 

contains the contact interaction potential ^(x, x') ~ 
-g<5(x-x'). 

We choose the disorder potential Vdis(x) to describe 
the presence of random impurities and to be indepen- 
dent of the hyperfine state, a choice that can be easily 
relaxed. As discussed in the introduction this choice is 
quite different from the cases of disorder introduced by 
laser speckles or incommensurate lattice potentials. In 
the present case, we consider impurity atoms to be heav- 
ier than the other type of atoms in the superfluid state, 
and it is also more useful to have the impurity atoms 
being fermions to prevent any possibility of low density 
superfluid behavior, as would happen for bosons in op- 
tical lattices at low filling factors. So for instance, one 
could choose to be the impurity atoms randomly dis- 
tributed in a three dimensional lattice that overlaps with 
a cloud of ^Li atoms, such that ^Li and ''"K only inter- 
act with each other without feeling the presence of the 
other atom's trapping potential. This means that the 
lighter atoms (e.g. ^Li) do not feel the presence of the 
optical lattice trapping the heavier atoms (e.g. ''°K), and 
conversely the heavier atoms do not feel the presence of 
the harmonic potential confining the lighter atoms, but 
the lighter and heavier atoms interact with each other. 
In this case, the interaction between the heavier impu- 
rity atoms and the lighter atoms can be expressed by an 
effective impurity potential 

Vdis{^)=J2vdF{^-r,), (3) 

i 

where rj are the locations of the impurities, and F{IV) — 
7r~^/^£^^ exp(— _R^/€^) is a Gaussian of width £d describ- 
ing the range of the impurity potential. Notice that Vd 
is a measure of the interaction between heavier impurity 
atoms and lighter atoms, and F(x — r^) is a measure of 
the square of the Wannier wavefunction |vj(x — ri)p of 
the heavier impurity atom. If the interactions between 
every heavier impurity atom and every lighter atom are 
the same and always repulsive, the effective impurity po- 
tential is smooth and strictly repulsive. 

The average over impurity positions leads to the cor- 
relator 

(Kiis(x)ydis(x')) = nKi^ - x'), (4) 
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where k = riiV^ with Ui being the impurity concentration, 
and the Kernel 



K{-K - x') 



1 



(27r)3/2^3 



exp 



,'^2 



(5) 



By inspection, only in the limit of — > 0, the correlator 
becomes (Vdis(x)Vdis(x')) = k(5(x — x'), reflecting white- 
noise correlations. 

Given that the Hamiltonian is fully described, we dis- 
cuss next the effective action of the system, and the cor- 
responding approximations to obtain the critical temper- 
ature of superfluid fermions from the BCS to the EEC 
regime in the presence of disorder. 



III. EFFECTIVE ACTION IN THE PRESENCE 
OF DISORDER 

The effective action during the evolution from BCS to 
BEC superfluidity and the calculation of the critical tem- 
perature in the absence of disorder were obtained long 
ago [l^ . and are briefly reviewed here to provide back- 
ground for the derivation of the effective action in the 
presence of disorder. 

We begin our discussion by recalling that the effective 
action in the absence of disorder can be derived by intro- 
ducing the pairing field A(a;) = ('0^(a:)V't(3^)) as a ther- 
mal average, where x — (x, r), with r being the imag- 
inary time. Upon integration of the residual fermionic 
degrees of freedom, the resulting action up to Gaussian 
order in the pairing field [l^ has two parts 



ScfL — So + ^Ga 



(6) 



The first term describes the action of unbound fermions 



So 



-C(k)/T 



(7) 



with ^(k) = e(k) — /i, where e(k) = k"^ /2m is the kinetic 
energy of a fermion of momentum k. The second term is 
the Gaussian action 



(8) 



due to the pairing field A((7) and its complex conjugate 
A(g) in momentum space, the former of which is just the 
Fourier transform of A (a;). Here, 



9 ^ 2(ig£ 



(9) 



is the fluctuation propagator, and q — (q, \q() is the four- 
momentum with g£ = 2'iTtT being the bosonic Matsub- 
ara frequencies. Here, the function Xi = tanh[^i/2r] 
describes the occupation of fermions with energy ^1 — 
^(k — q/2), and the function X2 ~ tanh [^2/2T] describes 
the occupation of fermions with energy ^2 = -I- q/2). 



To derive the effective action for a fixed configura- 
tion of disorder, we define the local chemical potential 
/i(x) = fJ. — Vdis(x). This corresponds to the local den- 
sity approximation and treats the effects of the disorder 
potential semiclassically. In the case of the impurity po- 
tential originated by heavier atoms randomly distributed 
as described above, this is expected to be a reasonable 
approximation, as the impurity potential is purely re- 
pulsive and non-confining such that bound states of the 
impurity potential do not exist. The situation for speckle 
potentials or other random potentials that possess bound 
states is quite different, and will be briefiy discussed later. 

Within the local density approximation, the effective 
action in the presence of impurity disorder becomes 



SoffiVdis] = 'S'o[Vdis] + -Spair [Vdis] , 



(10) 



where S'o[Vdis] is the action of unbound fermions in the 
presence of disorder given by 



dx 



,-«(k,x)/T 



(11) 



with ^(k, x) = e(k) — A'(x). The second contribution to 
S'eft has two terms, that is Spair — Sq^u + S4. The first 
term corresponds to Gaussian pairing fluctuations 



/dx 
— ^A(g,x)r-i(<z,Kiis)A(g,x), 



where A(g,x) is the pairing field, and 



(12) 



V 



X1+X2 



^ ^"'^^ 9 ' ^ W^i - a - 6 - 2Kiis(x))] 

(13) 

is the pair correlation function in the presence of disorder. 
Here, the atom-atom interaction g can be expressed in 
terms of the scattering length ag as 



V 
9 



mV 
47rac 



E- 

^ 2ek 



(14) 



In addition, the notation q — (q, ig^) represents the four- 
momentum, the function Xi — tanh [(^ -t- Vdis(x)) /2T] 
describes the occupation of fermions with energy 
= ^k — q/2), and the function X2 — 
tanh [(^ -I- Vdis(x)) /2T] describes the occupation of 
fermions with energy ^ = ^(k -|- q/2). 

The interaction between the pairing fields is described 
by the action 



S,^{TV)-' J dxdr^|A(x,T)|^ 



(15) 



where we used the definition A(x, t) = A(q = 0,x, r). 
Here, the coefficient 



^-E 



Y 



X 



(16) 
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represents the strength of the effective interaction be- 
tween pairing fields, and the notations X = tanh(^k/2r) 
and Y = sech^(^k/2r) have been used. 

An expansion of the pair correlation function 
r^^((7, Vdis) in the limit of low energy (small ig^, with 
the analytic continuation ig^ — > + i5), long- wavelength 
(small |q|), and small disorder potential leads to 

r-^(g, Fdi,) =a + + eydis(x) + dio. (17) 

The coefficients of the Taylor expansion are 



a{p,T) 



mV 
47ras 



E 



1 X 
2^ " 2^ 



(18) 



describing the order parameter equation in the absence 
of disorder when a{^,T) = 0, 



c(/i,r) = 



k 



X 



Y 



XY kl 



(19) 



describing the kinetic energy of the pairing field, 
d[n,T) = dR + id/, with d^ = Ek^/^Ck ^^'^ ^.i = 
\TrN{ep)^/ {STy/eF)\ <d{fi) describing the frequency de- 
pendence, and e{^i,T) = ~da/d^\T describing the 
strength of effective disorder potential. 

It is very convenient to rewrite the effective action for 
the pairing field S'pair into a more familiar form. Thus, 
next, we put Eq. PT)) back into the action of Eq. (IT^ . 
set |qp — )• — V^, A((7,x) — > A(q = 0,i(7£,x), and scale 
the pairing field A/\/dR — !• ^ such that the coefhcient of 
the real part of duj becomes one. A subsequent Fourier 
transformation into imaginary time, and the inclusion of 
the fourth order term in 5* leads to the pairing field action 

^pair [14is] = ^Gau [^^dis] + ^4 [Kli.] , (20) 

where ^cauiVdis] is defined in Eq. ([T^ and 

^4[Vdis] = {TV)-' J dxdr^|*(x,r)|4, (21) 

where 5, = b/dj^ is the effective interaction. In order 
to explicit out dissipative and non-dissipative temporal 
fluctuations, the pairing action is rewritten as 



Vir[^,*,^dis] = (TV)-' / dxdr/: 



where 



(22) 



Cnd + is the Lagrangian density that 



contains two contributions, 
dissipative part 



The first term is a non- 



Cnd{x,t) = ^ 



dr- 



2m, 



- Ai* + 7^dis(x) 



(23) 

and describes a generalized Gross-Pitaeviskii Lagrangian 
for the scaled pairing field ^' = \I'(x, r). The term con- 
taining the effective mass m* = dum/c is the kinetic 



energy of the pairing field, fi^ = —a/dn plays the role of 
the pairing field chemical potential, and 7Vdis(x) is the 
effective disorder potential, with 7 — e/dR. The second 
term is the dissipative contribution that comes from the 
term containing idjuj and reflects the decay of fermion 
pairs into unbound fermions (Landau damping), which 
in imaginary time takes the Caldeira-Leggett form 



Cd{x,t) 



2^ 



dr 



,|vI/(x,T)-vI/(x,r')| 



'M2 



(r - T'y 



(24) 



where A ~ di/d^. In our present discussion A 7^ for 
/i < and A = for > 0, since di is non-zero for 
positive /i, but vanishes for negative ^. This indicates 
that on the BCS side (/i > 0) the lifetime of pairs is 
short (pairs break) near the critical temperature, while 
on the BEG side /i < the lifetime of the pairs is infinite, 
or better said, very large, such that stable fermion pairs 
exist even above the critical temperature. 

It is important to emphasize that the long-wavelength, 
low-frequency effective action derived in Eq. (1^^ is valid 
for the entire evolution from BCS to BEC superfluid- 
ity near the critical temperature Tj., so long as the lo- 
cal density approximation is applicable. Now that we 
have derived the effective action for the pairing field in 
the presence of disorder, i.e., a purely repulsive random 
impurity potential, we discuss next the self-consistency 
equations that output the critical temperature as a func- 
tion of disorder for a given interaction parameter. 



IV. THE CRITICAL TEMPERATURE IN THE 
PRESENCE OF DISORDER 

In order to determine the critical temperature in the 
presence of disorder, we need to derive the correspond- 
ing order parameter and number equations. For this 
purpose, a natural choice of dimensionless parameters 
are l/lkpcis) for interactions, rj = Knp/ejp for disorder, 
and T = T/ep for temperature, where kp is the Fermi 
momentum, np = fc|,/37r2 is the fermion density and 
ep = kp/2m is the Fermi energy. Notice that our di- 
mensionless parameter rj that characterizes the degree of 
disorder is directly obtained from the original Hamilto- 
nian density in Eq. ([1]) and from the disorder correlator 
(Vdis(x)Vdis(x')), and is thus an input parameter of the 
theory. Since k = riiV^, and Vd has dimensions of energy 
times volume, it can be written as the product of the 
impurity volume -^^ and the characteristic amplitude Vd 
of the disorder potential Vdis(x), i.e. Vd = Vdf^d- There- 
fore, we may write t] = {ni£^){np£^){vd/eF)'^ , which re- 
veals that r] is small when the range of the impurity 
potential id is much smaller than the average separa- 
tion between impurities (n^ ^^'^ id ^ 1) or between 

fermions n^^^ {np^^^id ^ 1), and also when the ampli- 
tude of the disorder potential Vd is much smaller than the 
Fermi energy ep {vd/sp ^ !)■ In some of the older litera- 
ture discussing the effects of disorder for non-interacting 
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fermions, the dimensionless parameter is typically cho- 
sen to be l/{kp£), where i is the mean free path. Note 
that is a derived quantity, and not an input parameter of 
the theory. Even though the mean free path is a useful 
concept, when the Fermi system is strongly interacting it 
becomes increasingly more difhcult to define it, and also 
to calculate it. Thus we prefer to use our input param- 
eter 77 to describe the degree of disorder of the system. 
Later in our discussion, we will return to the connection 
between our dimensionless parameter 77, and the mean 
free path of paired and unpaired fermions. 

The order parameter equation can be immediately read 
from Eq. (|23p and is given by the condition 



fJ,*{T, l/{kFas)) = 



(25) 



corresponding precisely to the order parameter equation 
in the absence of disorder a{fi, T) — given in Eq. ([T8| . 
indicating that this equation is not explicitly affecte d by 
weak disorder, as required by Anderson's theorem [l5j. 
However, the determination of the critical temperature 
Tc{rf) and the chemical potential /i(77) as a function of 
dimensionless disorder 77 for fixed scattering parameter 
iKkpas) requires the simultaneous solutions of the order 
parameter Eq. ()25|) and the number equation, which is 
discussed next. 

To determine the number equation in the presence 
of disorder, it is necessary to average the thermody- 
namic potential r2(Vdis) rather than the partition func- 
tion fl7l|. Thus, to compute the disorder average (...) = 
/ W P\V]{. . .) , where the probability measure is 



P[V] = exp 



i y"dxdx'ydis(x)/^-i(x - x')Fdi.(x') 



(26) 

we just need the disorder correlator if (x — x'). The 
thermodynamic potential for a fixed configuration is 
f^(Vdis) = — rinZ(Vdis), and Z(Vdis) can be expressed 
as the product of two contributions 



Z{Vdis) — ^o(Kiis) X Zpair(Vdis) 



(27) 



The first term is the partition function for unbound 
fermions 



Zo(Kiis) = exp[-5o(Kiis)] 



(28) 



and the second term is the partition function for paired 
fermions 

^pair(14is) = exp [- Vh^^-, ^, Fdi.)] . (29) 

Thus, fi(Vdis) = f^o(Vdis) + f^pair(Vdis), whcrc the un- 
bound fermion thermodynamic potential is f2o(Vdis) = 
—Tin Zf){Vdia) = T SoiVdis) ^ while the contribution due to 
paired fermions is f^pair(Vdis) = -T'ln 2'paii.(Vdis)- How- 
ever, the calculation of the average (51) = -~T(ln Z) can 
be performed by using the standard replica trick [fTll for 
statistical averages, where 



(InZ) =limM^oln(Z*0'/*'. 



(30) 



This trick was first introduced in the context of spin 
glasses JJi] , and can be applied to both the unbound and 
paired fermion parts of the effective action. Here, how- 
ever, we discuss explicitly the application of the replica 
trick only to the more complex case of the pairing field, 
and just quote the result for the thermodynamic poten- 
tial of unbound fermions. In the case of the pairing field, 
the replicated (M-copies) partition function is 

„ A/ f M ^ 

j=l I 1=1 J 

(31) 

with S'pair defined in Eq. (1^^ . Taking the configu- 
rational average (.^pfij.) with the probability measure 
P [Vdis] amounts solely to a Gaussian integral over Vdis, 
which can be easily performed leading to 



7M 



V 



>i, exp ■ 



M 



(*,0) 



M 

E 



(32) 

Here, we used the notation I?^, = Ili^i ^ with 
S'pair («, 0) = 5pair [^i, ^i, Vdis = O] Corresponding to the 
pair action without disorder and with the replacement 
of * ^- *i in Eq. The second term in Eq. ((5^ 

corresponds to a density-density interaction between the 
replicated fields given by 



5dis 



Tdrdxdx pi{x)K — x')pj{x') 



(33) 

where Kpair = 7^k. The densities Pi{x) — 
^i(x, T)vl'j(x, r) and Pj{x') = ^'j(x', t)5'j(x', r) appear 
for equal time r, while K{x — x') is the disorder correla- 
tor defined previously. Thus, the replica trick transforms 
the problem of the pairing field in a disorder potential, 
into the problem of interacting replicated pairing fields 
^I^i with interaction strength — (Kpair/2)-R'(x— x'), in ad- 
dition to the standard direct interaction term character- 
ized by the interaction {g* /2) S{x — x'). There are now 
two types of interactions between the replicated pairing 
fields as indicated in Fig. [2j 

The calculation of the averaged thermodynamic poten- 
tial for paired fermions 



pair/ = ~T lim ln{Z 



M \1/M 



(34) 



can now be performed perturbatively as a function of 
the disorder parameter Kpair and the interaction g^. For 
this purpose it is better to work in four-momentum space 
q = (q, ig^) and obtain the replica Green's function 



GJiM) = -i<7£ + £^pair(g), 



(35) 



where impair (g) 



Mul + |qP/2m, - ^, -I- I]pair(<7)- Here, 
5^pair(<z) is the Dyson's self-energy for paired fermions 
containing all possible Feynman diagrams to the desired 
order in powers of 5* and /tpair- These types of diagrams 





FIG. 2: The two types of interaction lines between the repU- 
cated paired fields. The diagram in a) is the direct intra- 
replica interaction, containing g-t, while the diagram in b) is 
the inter-replica interaction containing k. 



are standard in diagrammatic perturbation theory [lOj. 
For instance, to first order in 5, and Kpair the self-energy 
has only 2x2 = 4 Feynman diagrams as shown in Fig. |3l 
while to second order in and Kpair there are 8 x 4 = 32 
diagrams as seen in Fig. S) 






P 



P 




a 



a 



FIG. 3: Two possible first-order self-energy diagrams side by 
side for the each interaction line of type a. Here, a describes 
either the intra-replica or inter-replica interaction line. 

The self-energy can be written as the expansion 



-'pair 



.(1) 
■^pair 



pair ' 



where Sp^j^. corresponds to 

(2) 



the diagrams in Fig. [3l and Spair diagrams in 

Fig. m The replicated partition function is {Z^ 

n,Gp^ir(9)/T 



-M 



M \ 
pair/ 



leading to the pair thermodynamic 



(-,-1 

pair 



(9)/T 



from which 



potential (repair) = T'Z^qln 
the contribution to the number of paired fermions can 
be obtained via the thermodynamic identity A^pair = 

The disorder averaged thermodynamic potential is 
(fi) — {Hq) + (repair), where the first term (also cal- 
culated using the replica trick) corresponds to the 
thermodynamic potential of unbound fermions (r^o) = 
-TEk^^AGo\k)/T] and 



Eoik) 



(36) 



is the fermion Green's function for unbound fermions 
with four- momentum k = (k, ikm) in the presence of the 
disorder potential. Here, Eo{k) = k'^ /2m — n + So(fc), 
where So(fc) is the self-energy for the unbound fermions 
in the presence of the disorder potential. This self- 
energy can be written as the expansion Eq — kOT^q /dK + 



P 



P 



FIG. 4: Eight possible second order self-energy diagrams side 
by side for fixed interaction lines labeled a and /3. Here, a and 
l3 describe either an intra-replica or inter-replica interaction 
line. 



^K^d^T,o/dK^ + 0{K^). The number of unbound fermions 
is then given by Nq — ~d{no) /d^i, and the full number 
equation is = — where 



7V(k) =iVo(K) + 7Vpair(K) 



(37) 



contains the contributions No{k) = 
TJ2kGo{k)dEo{k)/diJ, for unpaired fermions, and 
Npi,n-iK) = -T^;, Gpair(9)9£^pair(g)/9/i for paired 
fermions. 

In our approximation, the solution to Eqs. (j25p 
and (j37p produces the critical temperature Tc as a 
function of both disorder 77 and interaction parameter 
l/^kpa-s), as shown in Figs. [S] and [51 Our phase diagram 
for superfluidity is obtained at finite T and to second- 
order in K (rj), in contrast to T = properties 1^,1^ (such 
as superfluid density or condensate fraction) obtained to 
Hnear order in k (77). Although the replica technique can 
be used to all orders in k (77), we perform numerical calcu- 
lations only to second order in k (77), because it becomes 
impractical to calculate higher order contributions to the 
self-energies. Therefore, our critical temperatures can be 
expressed as 1^(77) = f;(0) [l - arj - Ptj"^] , where 1^(0), 
a and (3 are functions of l/{kpas) only. In Fig.O we show 
Tc/ep versus 77 for three different values of the interac- 
tion parameter l/ikFCis) corresponding to the BCS side 
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FIG. 5: The critical temperature Tc/tp versus dimension- 
less disorder parameter rj for different values of (fcFfls)"^ = 
— 1.2, 0, 1.2, respectively. 



at l/(fci?as) = —1.2, to unitarity at l/ikpag) = 0, and 
to the BEC side at l/ikpas) = 1.2. While in Fig. El we 
show a three-dimensional plot of T^/ep as a function of 
l/{kpas) and rj. Two important points can be inferred 
from Figs. [5] and |6] describing the critical temperature. 
First, impurity disorder effects seem to be more detri- 
mental in the BEC regime, when compared to the BCS 
regime. Second, it is near unitarity ((fci?as)~^ ~ —0.3) 
that the superfluid is more robust to impurity disorder. 
We attribute this qualitative difference to the fact that in 
the BCS regime phase coherence can only be destroyed 
through simultaneous breaking of pairs, which does not 
occur since the impurity scattering is elastic, and impuri- 
ties can not provide the required energy break the pairs. 
On the other hand, in the BEC regime the pairing field 
corresponds to molecular bosons, and phase coherence 
can be more easily destroyed without the requirement of 
simultaneously breaking pairs. 

Even though numerical results are useful, analytical 
insight into the evolution from BCS to BEC superfluid- 
ity in the presence of disorder can be obtained through 
calculations of the critical temperature to first order in 



1/ {kpa,) 

FIG. 6: Critical temperature Tc/ep as a function of dimen- 
sionless disorder r] and interaction l/(fc_Fas) parameters. 

77, as discussed next. 



V. ANALYTICAL RESULTS FOR WEAK 
DISORDER 

Analytical results to first order in the dimensionless 
parameter 77 can be obtained in both the BCS and the 
BEC limits, where simplifications occur and calculations 
are easier, but yet non-trivial. 

In the BCS limit the interaction parameter is large and 
negative {\/{kpas) ^ —1), and the number equation is 
dominated by the contribution due to unbound fermions 
iVo 3> A^pair- Given that the order parameter equation 
described in Eq. ((25)) is unchanged, then any change in 
the critical temperature Tc has to come from a change in 
the chemical potential due to the presence of the impuri- 
ties. This change can be computed through the number 
equation defined in Eq. (j37p by including only the domi- 
nant contribution due to unbound fermions. 

The self-energy to linear order in k is Eq — 
kI/-i ^j^(ia; - ^k)"^ ~ /^^(/i) [Ac - i7rsgn(w)] , where 
S{ji) = (3/2) Ac = vr^''''Ay^/{l - y), and 
^d/^F = ^/{kpid)'^ takes into account the range id of 
the disorder correlator K{x — x'). Replacing Eq into 
the number equation, and performing the integrations 
over four-momentum leads to a fermion density n = 
N/V nn + Sn, where no = 2Y,kf ~ Cji^^^ and 
dn = irjCK^jJL with C = {A/i)V-^M{€p)ep^'^ . In 
the case of perfect particle-hole symmetry 5n can be ne- 
glected and the chemical potential is pinned to the Fermi 
energy {ji = ep) leading to unchanged given by 

Tc = ej.(eV^)8e-2g^p[_^/(2fcj.|a,|)]. (38) 

This refiects Anderson's theorem, as changes in Tc can 
occur only via the disorder dependent shift in the chem- 
ical potential. Relaxing the condition of perfect particle- 
hole symmetry leads io ji — ep[\ — D7f\, where D = 
2 Ac — 4(ed/e_F)^/^. The change in /x produces a corre- 
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spending change in Tc{ri) = T'c(O) [1 — arj] , where a = 
TTD/{4\kFas\). 

In the BEC regime {l/{kpas) ^ 1) the number equa- 
tion is dominated by iVpair ^ -^o given that all fermions 
are paired into molecular bosons. To linear order in k the 
pair self-energy is Epairlg) = .9* "-pair/2 -t- Kpair^(g) [13) 
where A{q) = A(qc) — ImJ"^ j (47r)-y/|^*| — \qi, with 
^(Qc) = Z]|q|<|q^| ["^*/(|ctP^)] ^-nd the renormalized 
chemical potential is /I* = -f- .9*npair/2 -I- KpairA(qc), 
with ripair = N-p^^jV. Noticing that dEpa.ixld^i —2 
in the BEC limit, and upon summation over Matsubara 
frequencies iqi in Eq. ([37]), we arrive at ripair — '^ub, 
with UB = C(3/2)(TO*Tc/27r)3/2 + Kp^i.Tcml / Att'^ when 
jj,* = 0. The solution of the number equation gives the 
critical temperature 



T,{k) = T,(0) [1 - Kpah.Te(0)m^/6^2„B] , 



(39) 



where TdO) = 2t: [nB/Ci3/2)f^^ /m^ is the Bose- 
Einstein condensation temperature for a gas of molec- 
ular bosons, which in terms of rj leads to Tc(r]) — 
fc{0) [1 - m]] . When the BEC limiting values Tc(0) 
0.218ei?, TO* — )■ 2to and np/2 are used, the coeffi- 

cient a — l2Tr^(Tc{0)/eF) « 25.8 is large. It is very im- 
portant to emphasize that if temporal fluctuations were 
neglected in the self-energy Spair((z) by setting iqi — 
then one would have come to the conclusion that the crit- 
ical temperature Tc for the paired fermions (molecular 
bosons) is essentially unaffected by disorder. However, 
our calculations show that disorder affects the phase co- 
herence of the molecular bosons via the incoherent part 
of Gpair('?) manifested in the branch cut -^/l/x^l — iqi of 
the self-energy Spair(9)- It is these quantum (temporal) 
phase fluctuations, which lead to a strong suppression of 
Tc in the BEC limit (in comparison to the BCS regime), 
where fermions are largely non-degenerate, particle-hole 
symmetry is absent, and Anderson's theorem is not ap- 
plicable. 

Now that we have established analytically and numer- 
ically the critical temperature of a disordered superfluid 
from the BCS to the BEC limit as a function of the di- 
mensionless impurity parameter 77, which measures the 
disorder correlation energy with respect to the Fermi en- 
ergy, we would like to comment briefly on the relation 
between rj and the unbound fermion mean free path £p 
and the relation between rj and the pair mean free path 

^pair- 



VI. MEAN FREE PATHS 

The concept of mean free path £ is often used in con- 
nection with the classical idea of the average distance 
travelled by a particle between collisions with impuri- 
ties. For quantum particles, however it is necessary to 
formulate a more precise meaning for the mean free path 
in terms of the details of the impurity potential. For 
instance, consider the transition amplitude U(x.,y,t) — 



(y| exp(— ii/i)|x) for a quantum particle to propagate 
from position x to position y for a given realization of 
the impurity potential contained in the Hamiltonian op- 
erator H . This amplitude can be thought as the sum of 
all Feynman paths connecting the two positions, which 
implies that the action for each path depends strongly 
on the each impurity configuration. Taking the impurity 
average of the transition amplitude (J7(x, y, t))dis leads 
to the averaging of random scattering phases, and the ex- 
pectation that translational invariance is restored by the 
averaging process, followed by a rapidly decaying ampli- 
tude: (C/(x, y,f))dis exp(|x — y|/£). The decay con- 
stant £ is called the elastic mean free path. Notice that 
the very definition of £ requires an exponential decay of 
the (/7(x,y,i))dis. In the cases where this is not true, it 
is more difficult to define the elastic mean free path and 
thus to use it as a measure of the strength of the im- 
purity potential. It is important to emphasize that the 
concept is only useful if indeed one can show that this 
exponential behavior exists. 

In two simple limits, we can relate our dimensionless 
disorder parameter 77 to an appropriately defined elas- 
tic mean free path £. For this purpose, we look at the 
imaginary time description of the transition probability. 

In the BCS limit, the Fourier transformation of the un- 
bound fermion Green's function Go{k) defined in Eq. ([551) 
into real space and imaginary time leads to Go(x, y, r) = 
Go,no(x, y, r) exp [— |x — y|/^i^-], where £f is the elastic 
mean-free-path of unbound fermions, and Gno is the sin- 
gle particle Green's function in the absence of disor- 
der. In this case, £f can be calculated analytically as 
kp^F = 4/(37rry) to leading order in rj. This result shows 
that when rj is small then kF£F is large. 

In the BEC limit, it becomes more appropriate to 
look at the paired fermion mean free path. In this 
case, we look at the Fourier transform of Gpair('?) 
defined in Eq. ([55)) . which becomes Gpair(x, y,T) — 
Gpair,no(x,y,T)exp [-|x- y|/£pair], whcrc £pair is gen- 
erally a function of r, and Gpair.no is the pair Green's 
function in the absence of disorder. When only low 
frequency contributions are included we have ^pair = 



47r/ 



f-pair 



/2 



1/2 



which in units of k. 



^pair 

becomes kF£pair = F{l/kFas)/ ^Jfj. Here, F is just a 
function of l/kpas, which takes the limiting value of 
F{l/kFas 00) ~ 0.26. Notice the non-analytic be- 
havior of /cF^pair whcn 77 — > 0. 

In the BCS and BEC cases the phase boundary be- 
tween the superfluid state and the normal state occurs 
prior to the limit where kplp — 1, and fc^^pair = 1, re- 
spectively. However, the present estimates do not include 
the effects of Anderson localization, which may become 
important at low-temperatures {T 0) and at larger 
values of 77. 

Having made the connection between the dimension- 
less disorder parameter 77 and elastic mean free paths for 
repulsive impurity potentials, we would like to make next 
some comparative remarks between speckle and impurity 
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potentials. 

VII. DIFFERENCES BETWEEN SPECKLE AND 
IMPURITY POTENTIALS 

Before concluding, we would like to point out some 
key differences between speckle and impurity potentials 
to emphasize that not all disorders are equal. 

The impurity potential described here is purely repul- 
sive, meaning that it is non-confining, in such a way that 
for each configuration of disorder there are no confined 
single particle states. Thus, the impurity potential rep- 
resents essentially scattering centers which affect mostly 
the phases of the wave functions of the particles that are 
scattered from them. In particular in the BEC regime, 
phase-fluctuations of the pairing field occur more easily 
than in the BCS regime, and thus impurity effects are 
more detrimental relatively speaking. 

In the case of speckle potentials, the situation is quite 
different, as there are many mountains and valleys in the 
disorder potential landscape, which allow for confinement 
and the existence of spatially confined states. While in 
the BCS regime these differences between speckle and 
impurity potential are less dramatic due to the robust- 
ness of Cooper pairs to phase fluctuations, the situation 
in the BEC regime is different. In the BEC limit the 
speckle disorder potential seen by bound pairs allows 
Bose-Einstein condensation into the spatially confined 
states creating many lakes (local condensates) near the 
minima of the disorder potential. As the average ampli- 
tude of the speckle potential increases lakes are formed 
and become disconnected beyond a certain threshold. 
Once the connections between the local condensates is 
lost, global phase coherence through the entire cloud is 
lost, and superfluidity disappears leading to an insulat- 
ing state. Thus, the problem of disorder due to speckle 
potentials is more like a percolation problem, and has 
a percolation threshold. While the confining regions of 
the speckle potential favor Bose-Einstein condensation, 
overlap between regions of local condensation determine 
the global phase coherence and thus superfluidity in the 
system. Because of this qualitative difference that fa- 
vors local Bose-Einstein condensation, which is a neces- 
sary but not sufficient condition for Bose superfiuidity 
in three-dimensions, we expect that speckle potentials 



have less impact on the critical temperature for super- 
fiuidity (as found in quantum Monte Carlo simulations 
for atomic Bose systems in speckle potentials than 
the purely repulsive impurity potentials described here. 
A more quantitative comparison between the effects of 
disorder due to speckle potentials and repulsive impurity 
potentials requires modifications of the current approach 
in order to include the very important contributions of 
confining regions, which are of course absent in the repul- 
sive impurity potential treated here. In addition, one has 
to be careful in comparing the sizes of molecular bosons 
(produced via Feshbach resonances) and the character- 
istic lengths of speckle potentials to determine what is 
the effective disorder potential felt by such molecules. 
These important theoretical comparisons between differ- 
ent types of disorder are left for a future publication, 
but we urge the experimental community to explore the 
differences between these types of potentials. 



VIII. CONCLUSIONS 

In conclusion, we analyzed the effects of repulsive im- 
purity disorder potentials on the critical temperature for 
superfluidity of ultracold fermions during the evolution 
from the BCS to the BEC regime. For s-wave superfluids, 
we showed that weak disorder does not affect the critical 
temperature of a BCS superfluid with perfect particle- 
hole symmetry in accordance with Anderson's theorem, 
as the breaking of fermion pairs and the loss of phase co- 
herence occur at the same temperature. However in the 
BEC regime, phase coherence is more easily destroyed by 
repulsive impurity disorder without the need of simulta- 
neously breaking fermion pairs. Thus, in the BEC regime 
a more dramatic change in the critical temperature oc- 
curs, when compared to the BCS limit. Finally, we also 
showed that the superfluid is more robust to the presence 
of disorder in the intermediate region near unitarity. 
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